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Formation or destruction of hyperbolic chaotic attractor under parameter variation is considered 
with an example represented by Smale-Williams solenoid in stroboscopic Poincare map of two 
alternately excited non- autonomous van der Pol oscillators. The transition occupies a narrow but 
finite parameter interval and progresses in such way that periodic orbits constituting a " skeleton" of 
the attractor undergo saddle-node bifurcation events involving partner orbits from the attractor and 
from a non- attracting invariant set, which forms together with its stable manifold a basin boundary 
of the attractor. 
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I. INTRODUCTION 

Uniformly hyperbolic strange attractors have strong 
chaotic properties and are structurally stable, i.e. in- 
sensitive to variations of functions and parameters in 
the dynamical equations [H-Q- Formal examples are 
Smale-Williams solenoid, Plykin attractor, and some 
other mathematical constructions suggested mainly in 
1960th-1970th [3|-l6|. Recently a number of physically 
realizable systems with hyperbolic attractors were pro- 
posed [1CH13I]. 

One promising direction of search for real-world non- 
linear dissipative systems with the hyperbolic chaotic at- 
tractors may be based on consideration of scenarios of 
their appearance under variation of control parameters. 
(Note that this is an issue certainly different from the 
commonly known scenarios of transition to chaos, like 
that of Feigenbaum, which lead to non- hyperbolic attrac- 
tors.) Possible appearance of hyperbolic chaotic attrac- 
tors was discussed by Ruelle and Takens in the context 
of the onset of turbulence [H], EH , but they did not con- 
sider concrete examples of dynamical equations demon- 
strating the phenomenon. More recently, Shilnikov and 
Turaev [16|, E3 have suggested a scenario of emergence 
of attractor of Smale-Williams type in a kind of the so- 
called blue sky catastrophe. An example of dynamical 
equations manifesting this scenario was presented and 
studied numerically in Ref. [18[ . 

In this paper we discuss an alternative scenario for the 
onset or destruction of attractor of Smale-Williams type 
associated with collision of two chaotic invariant sets, an 
attracting and a non-attracting one. Concretely, we con- 
sider the phenomenon in a physically realizable system 
composed of two alternately activated non-autonomous 
self-oscillators, which pass the excitation each other in 
such way that the phases of oscillations at successive 
stages of activity evolve chaotically in accordance with 
an expanding circle map (Bernoulli map) [MEH- In this 
system attractor of Smale-Williams type occurs in the 



stroboscopic Poincare map. At certain parameters the 
uniformly hyperbolic nature of this attractor was verified 
in computations [HIEilES? 20]. Here we consider birth or 
disappearance of the chaotic attractor under variation of 
a parameter controlling a relative duration of the stages 
of activity and silence. The scenario may be thought of as 
a multitude of pairwise collisions of orbits relating to the 
attractor with those from some non-attracting invariant 
set happening in a narrow but finite parameter interval. 
In Appendix we shortly discuss a transparent example 
of a simple model map demonstrating qualitatively the 
same scenario of birth and destruction of chaotic attrac- 
tor. 



II. THE BASIC MODEL 

The system we will study is governed by a set of dif- 
ferential equations EH 



UJqU, 



x 2 )u — uoqx + — ycosujot. 



u = (h + Acos^ 
y = 2uj v, 

v = (h-Acos^- y 2 )v - 2uj y + ^x 2 , 



(1) 



slightly modified by introducing an additional parame- 
ter h. The model is composed of two subsystems, the 
van der Pol oscillators with characteristic frequencies ujo 
and 2cjo- The dynamical variables x and y are the gen- 
eralized coordinates of the oscillators while u and v are 
the normalized velocities. In each oscillator the parame- 
ter responsible for the birth of the limit cycle is forced to 
vary slowly with period T and amplitude A, in opposite 
phases for the two subsystems, which become active turn 
by turn. The newly introduced parameter h controls the 
relative duration of the stages of activity and damping. 
The coupling between the subsystems characterized by 
coefficient e is established in such special way that the 
excitation transfer between the subsystems is accompa- 
nied by doubling of the phase shift attributed to the os- 
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cillations on each next cycle of the transfer. Due to this, 
the stroboscopic map of the system (pQ) in a wide param- 
eter range gives rise to a hyperbolic chaotic attractor of 
S male- Williams type as it was verified in computations 
based on the cone criterion [12|,[l9| and on the statistics of 
angles between the stable and unstable manifolds 
At a particular parameter set ft = 0, ojq = 27r, T = 6, 
A = 5, e = 0.5 the uniformly hyperbolic nature of the 
attractor was confirmed with a mathematically rigorous 
computer-assisted proof [20j . 



III. NUMERICAL RESULTS AND DISCUSSION 




Let us start from the case ft = 0, where the opera- 
tion mode corresponding to the Smale- Williams attrac- 
tor certainly takes place, and move to the domain ft < 
increasing the absolute value of ft. Then, duration of the 
activity stages will decrease, and at some place the trans- 
ferred excitation level becomes insufficient to recover the 
amplitude of the oscillations. Then, the system drops af- 
ter some transient to the trivial stable state of zero am- 
plitude. Figure Q] shows the Lyapunov exponents of the 
Poincare map versus parameter ft. The right-hand part 
of the plot corresponds to situation of existence of chaotic 
uniformly hyperbolic attractor of Smale- Williams type. 
A characteristic feature is that the largest Lyapunov ex- 
ponent is positive, close to In 2 « 0.693 while others are 
negative, smoothly depending on the parameter. In fact, 
at ft < this attractor coexists with another one, the 
trivial stable zero state. (It is so because the coupling 
term in the last equation (pQ) has quadratic dependence on 
the dynamical variable; thus, with initial amplitude small 
enough, the system inevitably relaxes to zero state.) As 
seen from Fig.[TJ with decrease of ft, the breakdown of the 
chaotic attractor occurs at certain parameter value. As 
the remaining attractor is the fixed point in the origin, 
after the breakup of the chaotic attractor all trajectories 
tend there. 

The sudden destruction of the chaotic attractor indi- 
cates the crises-like (" dangerous" ) character of the transi- 
tion and makes one to suspect that the event is associated 
with a collision of two invariant sets (cf. |2l|). To bring a 
qualitative basis for this observation let us develop some 
approximate approach to stroboscopic description of the 
dynamics. It is convenient to introduce polar coordinates 
on the plane of the dynamical variables of the first oscil- 
lator. We set {xo,^i} = {x,u/0.9 — x/2} and introduce 
the amplitude (or radial component) and the phase (or 
angular coordinate) in such way that xq + ix\ = re w . 
(In definition of the variable change the coefficients are 
chosen to get the form of the attractor close to a circu- 
lar one in the plane projection.) Now, we construct a 
one-dimensional amplitude map for the first oscillator as 
follows. Starting with zero amplitude of the second oscil- 
lator and assigning some amplitude r and phase ip to the 
first oscillator, we evaluate the amplitude r new after one 
period of the parameter modulation and plot it as a func- 



FIG. 1: The Lyapunov exponents of the system (pQ) versus 
parameter ft (black curves) and the Lyapunov exponent for 
the trivial attractor, the fixed point at the origin (gray line). 
Parameter values are cjo = 27T, T = 6, £ = 0.5. 




FIG. 2: Constructing the amplitude stroboscopic map at 
loo = 27T, A = 6.5, T = 6, e = 0.5 for different values of 
parameter ft: (a) the original diagram accounting all accumu- 
lated computational data, and (b) the idealized plot excluding 
the widening by the phase averaging. 



tion of the initial amplitude. The constructed map is not 
so bad tool for sketchy description of the stroboscopic 
amplitude dynamics: since the subsystems are excited 
alternately, each epoch of activity for one of them cor- 
responds to silence period for another one, whose ampli- 
tude in the rough approximation can be neglected there. 

Figure [2] shows the plots for the amplitude map at four 
different values of ft as obtained from numerical integra- 
tion of equations (pQ). Actually, instead of curves we get 
the widened strips on the panel (a); the reason is that 
the positions of the plotted points depend not only on 
amplitudes but also to some extent on the initial phases, 
varied from to 2tt. Roughly, this widening measures 
a degree of incorrectness of the description in terms of 
the amplitude map. On the panel (b) this widening is 
excluded by means of averaging over the initial phases. 

As seen from Fig. (2}d, with variation of ft a tangent 
bifurcation occurs. At larger ft the map possesses three 
fixed points (on the plot they correspond to crossings of 
the curve with the bisector). One is located in the ori- 
gin that is the stable stationary zero state O; two others 
correspond to some finite amplitudes. The fixed point A 
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FIG. 3: Diagrams illustrating dynamics of phase on successive 
stages of excitation of the first oscillator on the attractor A 
(black) and on the non-attracting invariant set S (gray). 



at larger amplitude is stable while that with less ampli- 
tude S is unstable. As value of h decreases, the fixed 
points A and S approach each other, merge (approxi- 
mately at h « —1.23) and then disappear. Of course, 
this is only a preliminary rough picture, and we must 
discuss now what happens actually nearby this parame- 
ter value in the original system. 

The fixed point A for the original system (pp) is as- 
sociated with the Smale-Williams attractor: while the 
amplitude evolves almost periodically (with the period 
of the parameter modulation) the phases at successive 
periods behave in accordance with the chaotic Bernoulli 
map. Of the same kind is dynamics on a non- attractive 
invariant set corresponding to the unstable fixed point S; 
the only difference is the additional instability of orbits in 
this set in respect to the amplitude variations, as evident 
from Fig. [2)3. Henceforth, we use the designations A for 
the attractor and S for the non-attracting invariant set of 
the original system rather than for the fixed points of the 
amplitude map. Figure [3] illustrates the phase dynamics 
on these sets with the computed iteration diagrams. The 
branches shown in black relate to the attractor, and the 
gray ones to the non- attractive invariant set. For the at- 
tractor the plot is obtained in straightforward way in the 
course of numerical integration of the equations (pQ) fiol — 
H2j[H1: we evaluate the phases cp n = aig[xo(t n )-\-ixi(t n )] 
from the instant states for the first oscillator at t n = nT 
and plot ifn+i versus ip n . To do the same thing for 
the non-attractive invariant set S special measures are 
needed to keep the trajectory on it; actually, the picture 
was obtained from a collection of data for periodic orbits 
belonging to the set. 

Accounting that the dynamics for the phases on both 
invariant sets A and S correspond to topologically equiv- 
alent Bernoulli maps, the orbits belonging to each of 
them are in one-to-one correspondence with a set of infi- 
nite two-symbol (binary) sequences. 



It is well known that a chaotic attractor contains a 
dense set of unstable periodic orbits which constitutes 
a kind of "skeleton" for the attractor [22]. In terms of 
symbolic dynamics, it corresponds to a class of the orbits 
associated with periodic binary sequences. This "skele- 
ton" is a convenient object to comprehend the dynam- 
ics on the whole attractor A. The same is true for the 
non-attracting invariant set 5, which possesses its own 
"skeleton" of unstable periodic orbits. The both sets are 
in one-to-one correspondence (as follows from identical 
symbolic dynamics): each periodic orbit on the attractor 
has a partner ("dual orbit") belonging to the invariant 
set S. Dealing with the "skeletons" it is possible to look 
over all periodic orbits (at least up to some reasonably 
large period) and analyze their bifurcations in depen- 
dence on the parameter h. 

A periodic orbit of the stroboscopic map belonging to 
the attractor A has a three-dimensional stable manifold 
and a one-dimensional unstable manifold (actually, the 
last one coincides with the whole attractor). The Floquet 
multiplier responsible for the instability may be evalu- 
ated as fif « 2 P , where p is the period of the orbit. The 
next multiplier, associated with the amplitude perturba- 
tion \±2 is less than 1. (It may be estimated as 
where / is the function plotted in Fig. [JJd, and ta des- 
ignates the position of the stable fixed point.) On the 
other hand, a periodic orbit belonging to the invariant 
set S possesses a two-dimensional stable manifold and a 
two-dimensional unstable manifold; one multiplier indi- 
cating instability in respect to the phase is /if « 2 P , and 
another one, associated with the amplitude instability 
/if, is essentially smaller: 1 < /if < A*f • Under decrease 
of parameter /i, mutually dual orbits, one from the set A 
and another from the set S undergo a "fold" (or tangent, 
or saddle-node) bifurcation: they become closer, coalesce 
at some bifurcation value /ibif and then disappear. 

If the description in terms of the amplitude map of 
Fig. [2] be exact, the bifurcations would happen simul- 
taneously for all pairs of the dual orbits, but actually 
it is not so: the bifurcation values /^if occupy a fi- 
nite, although narrow interval. (Roughly, this interval 
— 1.22 < h < 1.24 may be thought as associated with 
a finite width of the "curves" visible in the panel (a) 
of Fig. El) 

The drawn picture is supported by computations with 
Eqs. (pQ): we have traced the bifurcations for a large num- 
ber of pairs of dual orbits associated with periodic sym- 
bolic codes. Table 1 collects the bifurcation values huf 
for the orbits of period from 1 to 5 along with their binary 
codes. 

According to the bifurcation theory [23[, the fold bi- 
furcation may be regarded as corresponding to a turning 
point of a curve, which represents a pair of periodic orbits 
in the extended space, which is the state space comple- 
mented with the coordinate axis of the parameter h. This 
interpretation reflects the fact that for h > /ibif each sad- 
dle cycle embedded in the attractor has, as a counterpart, 
a dual saddle cycle from the non- attracting invariant set. 



4 



TABLE I: Bifurcation values of parameter h for the unstable 
periodic orbits of the system (pQ) at cjo = 27T, A = 6.5, T = 6, 



p 


symbolic codes 


ftbif 


1 


R 


— 1.22583188 


2 


LR 


-1.23358210 


3 


LRR 


-1.22879290 




RLL 


-1.23489088 


4 


LRRR 


— 1.22695959 




LRRL 


-1.23226821 




RLLL 


-1.23371422 


5 


LRRLR 


-1.23070842 




RRRRL 


-1.22628724 




LRRRL 


-1.23012599 




RLLLL 


-1.23249362 




LRRLL 


-1.23263428 




RLLRL 


-1.23467205 



At the bifurcation point they both merge and annihilate, 
or, to put it another way, the respective branches of the 
curve continue each other. 

Figure [3^ plots coordinates of points of the periodic 
orbits of periods p < 10 for the stroboscopic map of the 
system (pQ) versus parameter h in the three-dimensional 
extended phase parameter space {xo, xi, h}. On this pic- 
ture attractor corresponds to the outer part of the tube- 
like formation, while the non-attracting invariant set cor- 
responds to the inner part. It looks as going inside out 
in order to be eventually transformed into the Smale- 
Williams attractor after passage of the turning points. 
Figure [3b plots radial (amplitude) component for the 
same set of periodic orbits versus parameter h; it may be 
regarded as a longitudinal slice of the object of panel (a). 

The structure observed at fixed h = —1.15 is illustrated 
in Fig. [5l Panel (a) represents, in fact, a transversal sec- 
tion of the object shown in Fig. [4^. Observe two forma- 
tions looking in rough approximation as closed circular 
curves but actually having fine transversal fractal struc- 
ture. The outer one represents the chaotic attractor A, 
while the inner one is the non-attracting invariant set S 
consisting of saddle orbits dual to those embedded in the 
attractor. It should be remembered, that the picture ac- 
tually is a projection of the corresponding sets from the 
four-dimensional state space of the stroboscopic Poincare 
map onto the plane. In Fig. [5]o the three-dimensional 
version of the diagram is shown, and one may see more 
details of the fractal-like transversal structure of both 
sets A and S. Evidently, the invariant set S and its sta- 
ble invariant manifold separate the basin of attraction 
of the Smale-Williams solenoid and that of the attrac- 
tive fixed point in the origin. (In Fig. [5^ the last one 
corresponds roughly to gray area inside the inner closed 
"curve" .) 

Objects similar to the invariant set S are encountered 
in the theory of complex maps being known as Julia 
sets [24], [HI- This analogy is rather deep: far from the 
bifurcation transition, when the size of the set S becomes 
small enough one can neglect all nonlinearities in the 
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FIG. 4: (a) Periodic solutions of the equations (pQ) with period 
< 7 plotted in the extended phase parameter space at cjo = 
27T, A = 6.5, T = 6, e — 0.5. (b) Superimposed graphs 
of radial (amplitude) component of the same set of solutions 
versus parameter h; this is actually a longitudinal slice of the 
figure (a). 

equations, beside that in the coupling term. Then, in 
a frame of description of the oscillators in terms of com- 
plex amplitudes, this set really turns to a Julia set for 
some complex map [261 ]. 

Qualitatively, the content of the considered scenario is 
illustrated by a simple artificially constructed model map 
in Appendix (in a similar way like one-dimensional maps 
illustrate the Feigenbaum or intermittent scenarios). 

IV. CONCLUSION 

Understanding scenarios for appearance of uniformly 
hyperbolic attractors under variation of control param- 
eters seems relevant for a search for real-world systems 
with such attractors, which will be of special practical 
and theoretical interest because of structural stability of 
the generated chaos. In this article, basing on numerical 
computations, we reveal a nature of the bifurcation sce- 
nario of birth or destruction of a uniformly hyperbolic at- 
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FIG. 5: Attractor A and non- attracting invariant set S in the 
stroboscopic Poincare section shown on the plane of variables 
of the first oscillator (a) and in three-dimensional plot (b) 
accounting a variable of the second oscillator at cjo = 27T, 
A — 6.5, T = 6, e = 0.5, and h — —1.15. The gray color on 
the panel (a) indicates roughly the basin of attraction of the 
stable fixed point in the origin. 



tractor of S male-Williams type in stroboscopic Poincare 
map of two alternately excited non- autonomous van der 
Pol oscillators under variation of a parameter controlling 
relative duration of stages of activity and silence of the 
oscillators. This is the second discussed in the literature 
scenario relating to formation of a uniformly hy per bolic 
attractor, after that of Shil'nikov and Turaev [16l-[l8j cor- 
responding to a kind of the so-called blue-sky catastro- 
phe. In our scenario, the birth or destruction of the 
Smale- Williams solenoid appears not as a single bifur- 
cation event, but occupies a finite parameter interval; it 
may be thought as a multitude of saddle-node bifurca- 
tions each involving a pair of orbits, one from the attrac- 



tor, and another from the non- attracting invariant set. 
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Appendix 

It is a good way to illustrate the considered scenario 
with a simple model. Let us start with the equation: 



Rzl 



(2) 



where z is a complex variable, and R is a real parameter. 
Setting z = re lkp we rewrite the map as 



fn+l 



= R- 



(f n +i = 2(p n (mod27r). (3) 



At small R the map for the radial variable r has only 
a single fixed point at zero. With increase of R at 
some instant (R = \/2), a tangent bifurcation occurs 
accompanied with appearance of a pair of fixed points 



\R? =p \ \R a — 1, one stable and another un- 



r l,2 = 

stable (see Fig. [6j). For the complex map (pp) the stable 
fixed point corresponds to attractor A placed on a cir- 
cle of radius r2 while the angular coordinate is governed 
by the Bernoulli map <^ n +i = 2cp n (mod27r). The un- 
stable point corresponds to an unstable invariant curve 
that is a circle of radius n on which the angular variable 
evolves in accordance with the same Bernoulli map, and 
this curve is a boundary separating basins of two coex- 
isting attractors, the stable point O and the attractor A. 
Now, consider a slightly modified version of the map 



RZnjZn 

^l + \z n (z n + s)\ 2 ' 



(4) 



which reduces to (0) at e = 0, and look what happens 
under increase of parameter R with e ^ 0. At small R 
the unique attractor is a zero fixed point O. For R great 
enough there is another attractor A, situated in the re- 
gion of relatively large \z\. Concerning dynamics on this 
attractor, the additional terms in the map (j4]) do not 
violate the nature of the dynamics of the angular vari- 
able, which follows a map of the same topological type 
as the Bernoulli map. However, now the invariant set A 
does not coincide precisely with the circle, but acquires 
a transversal split and looks similar to a projection of 
the Smale-Williams solenoid. The border between the 
basins of attractors O and A is represented now by some 
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r n 2 -2 Rez 2 



FIG. 6: A plot of radial component of the model map J2} (a) 
and disposition of attractors A, O and of a non- attractive 
invariant set S on the plane of complex variable z (b). 

non-trivial invariant set 5, which may be computed using 
backward iterations of the map (j4]): 

where a sign at the square root is selected on each next 
iteration randomly. Visually, the set S looks similar to a 
Julia set of complex quadratic map. (It is not surprising: 
neglecting the nonlinear term in the denominator of (j4]) 
we get just the complex analytic map z n +i = R(z^-\-ez n ) 
providing a reasonable approximation in the domain of 
location of the set S while \z\ is small.) Dynamics on 
the sets A and S may be described with the two-symbol 
Bernoulli scheme. (For the set S it is obvious from the 
above algorithm of backward iterations: a symbolic se- 
quence of an orbit corresponds to the sequence of signs 
selected in the course of the procedure.) 

With decrease of parameter R a size of the attractor A 



reduces while a size of the set S grows (Fig. [7j). As the 
sets A and S approach each other, the transverse struc- 
ture of S develops in the radial direction, and becomes 
less akin to the original Julia set. For some R = R\ a 
touch of A and S happens, and it corresponds to the mo- 
ment of crisis of the attractor A (the touch of the basin 
boundary, cf. |2l|). 

Each of the sets A and S contains a dense set of unsta- 
ble periodic orbits. A periodic orbit on the attractor A 
has a multiplier fif- = 2 P > 1, where p is a period of the 
orbit, and a multiplier p^ < 1. For an orbit on the in- 
variant set S the multipliers are /if = 2 P > 1 and /if > 1 • 
As symbolic representations of the orbits in A and S are 
of the same type, each orbit in A can be related to a 
partner dual orbit in S and vice versa. With decreas- 
ing R the dual orbits approach each other, then merge 
and disappear via the fold (tangent) bifurcation, where 
/i^ = /if = 1- For different periodic orbits this bifurca- 
tion occurs at different values of i?, so that the process of 
disappearance of orbits remaining on the set A occupies 
a finite parameter interval [R2,R\\. Only in the degen- 
erate case e = 0, all orbits belonging to the sets A and S 
undergo the bifurcations simultaneously, just at R = y/2. 

The above consideration is in clear qualitative corre- 
spondence with that in the main part of the article. Some 
inconsistence is that the model map (j4]) is non-invertible 
(in contrast to the Poincare map of the system (pQ)). This 
issue could be repaired with a Henon-like modification 
of the model: z n +i = Rz n (z n +e) _ jj Z D1 rt in 

this case the analysis and explanations become much less 
transparent. (We can keep up appearances saying that 
we deal with this invert ible map, but assign exclusively 
small value to the parameter b.) 
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